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ABSTRACT 


A rather  general  method  is  given  to  construct  classes  of  functions  with 
an  arbitrary  almost  constant  (linear)  initial  interval  followed  by  a non- 
prescribed  interval  of  monotonic  nonlinear  behavior.  This  region  of  non- 
linear behavior  is  succeeded  by  an  unbounded  interval  of  almost  constant 
(linear)  behavior.  They  contain  not  more  than  four  selectable  parameters 
and  are  synthesized  from  analytic,  monotonic,  normalized  and  bounded  base 
functions  through  the  introduction  of  two  separate  kernel  sets,  subsequent 
addition  and  integration.  As  examples  we  give  the  special  functions  based 
on  the  error,  the  hyperbolic  tangent,  the  inverse  tangent,  a rational  and  the 
incomplete  gamma  function.  Limiting  function  forms,  such  as  the  bilinear 
form,  are  derived  for  limiting  values  of  the  parameters . 


We  have  found  these  functions  useful  in  the  total  strain  approach  to 
viscoplasticity,  i.e.,  the  analytical  modelling  of  stress-strain  diagram.s, 
strain  (stress)-rate  effects,  creep  and  relaxation  curves  for  monotonic  and 
cyclic  loading.  Also  these  functions  offer  great  flexibility  in  the  curve 
fitting  of  experimental  data  generated  in  the  above-mentioned  tests. 


Introduction 


In  an  attempt  to  model  time  (rate) -dependent  as  well  as  time  (rate)- 
independent  behavior  of  metals  without  the  use  of  a yield  surface  and  without 
the  normal  decomposition  of  the  small  strain  tensor  into  elastic,  plastic  and 
creep  parts  a new  total  strain  approach  to  inelastic  behavior  was  proposed 
in  [1]  after  a careful  review  of  the  predictive  capability  of  currently  used 
constitutive  equations  vis  a vis  experimental  evidence  [2]. 

The  new  approach  employs  total  strain  only  and  recognizes  four  different 
phenomena  which  require  different  repositories  in  a constitutive  equation. 

They  are 

• nonlinearity 

• rate  (time) -dependence 

• history  dependence  in  the  sense  of  plasticity 

• aging. 

Rate  (time ) -independence  is  obtained  as  a special  case  of  rate  (time ) -dependence 
Aging  shall  be  excluded  from  our  considerations  in  the  present  paper. 

The  repository  for  the  modeling  of  history  dependence  in  the  sense  of 
plasticity  rests  with  a microstjructure  memory  function  which  can  discontinu- 
ously  change  at  points  of  unloading  [1].  Between  successive  points  of  unloading 
(as  defined  in  [1] ) the  constitutive  equ£.tion  is  a nonlinear  elastic  and  non- 
linear viscous  law  for  the  rate-independent  and  the  rate-dependent  case, 
respectively.  The  stress  and  strain  tensor  at  points  of  unloading  are  intro- 
duced into  the  constitutive  equation  to  affect  initial  elastic  response  upon 
unloading.  The  sign  of  the  differential  of  the  second  invariant  of  the  stress 
(strain)  tensor  measured  from  the  last  point  of  unloading  is  the  loading  and 


unloading  criterion  [1] . 


2. 

A total  strain  approach  to  plasticity  or  viscoplasticity  must  be  cJsle  to 
reproduce  the  almost  linear  behavior  of  metals  for  small  strains  over  an  almost 
arbitrary  range  followed  by  the  strongly  nonlinear  behavior  at  the  onset  of 
plastic  flow.  Moreover  since  the  "elastic  range"  may  be  changed  as  a function 
of  prior  deformation  (cold-working,  see  Pig. 4a  in  [3])  or  of  temperature  the 
functions  used  in  a total  strain  approach  must  be  capable  of  representing  these 
changes.  Our  approach  requires  that  this  be  done  by  merely  changing  the  con- 
stants in  the  functions  used  without  altering  the  function  form. 

An  example  of  such  a function  is  given  in  [3]  for  the  rate- independent 
case,  whereas  the  rate-dependent  case  is  treated  in  [4] . 

It  is  the  purpose  of  this  paper  to  give  ' rather  general  method  of  con- 
structing functions  of  almost  constant  (linear)  behavior  over  an  arbitrary 
interval  followed  by  nonlinear  monotonic  behavior  over  an  adjacent  interval. 
Finally,  almost  constant  (linear)  behavior  is  obtained  after  the  region  of 
nonlinearity.  We  propose  a whole  family  of  functions  in  closed  form  using 
not  more  than  four  selectable  constants.  Depending  on  the  value  of  these 
constants  the  value  of  the  ftinction  can  change  within  wide  limits  which  make 
these  functions  very  useful  in  their  applic.-tion  to  viscoplasticity. 

Construction  of  Nonlinear,  Differentirble  Functions  with  Selectable 
Intervals  of  Almost  Constant  or  Linear  Behavior 

An  examination  of  real  stress-strain  diagrams  and  the  required  properties 
of  functions  used  in  our  constitutive  theory  [1,3,4]  indicates  the  need  for 
the  development  of  a flexible  function  Y(X)  shown  together  with  Y'  (X)  in 
Fig.l  which  fulfills  the  following  requirements. 

The  function  Y(X)  must  be  odd  in  X,  with  Y' (0)  = S^;  Y(X)  must  be  linear 
in  X over  an  initial  interval  |x|  < X^,  followed  by  a non-prescribed  interval 
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3. 

of  monotonic  nonlinear  behavior.  Finally  an  unbounded  interval  |x|  > follows 
on  which  Y(X)  must  again  and  thereafter  be  linear  in  X with  constant  slope  S,. 

We  must  of  course  have  X^  > X^  > 0 and  while  we  are  primarily  interested  in 
the  case  >0,  ^ consider  any  other  finite  values  of  and  S^. 

We  choose  to  develop  a representation  of  Y(X)  which  is  a real  valued 
analytic  function  in  X for  all  finite  X . We  turn  to  the  development  of  Y- 
functions  which  are  approximately  linear  over  the  selectable  intervals  |x|  < 
and  |x|  > X^.  The  degree  of  approximation  is  controlled  by  a quality  factor 
which  can  be  arbitrarily  selected. 

We  accomplish  our  synthesis  of  Y(X)  indirectly,  by  first  constructing  Y' (X) 
and  then  integrating  in  X to  obtain  Y(X).  We  obtain  Y'  (X)  from  the  sum  of  a 
constant  plus  two  functions  such  that  on  specified  intervals  the  net  contribu- 
tion due  to  the  two  functions  is  approximately  zero,  leaving  only  the  constant 
behavior.  For  this  purpose  we  now  introduce  the  function  F(X);  let  F(X) 
satisfy  the  following  requirements  -. 

F(X)  must  be  a real  valued  analytic  function  of  the  real  variable  X 
P(X)  must  be  odd  in  X 

F(X)  must  be  bounded  at  X =*  co  and  normalized  so  tliat  F(oo)  =1 
F(X)  must  be  monotonic  increasing  as  X increases 

F' (X)  must  be  finite,  positive,  and  monotonic  decreasing  for  non- 
negative X. 

The  form  of  F(X)  is  that  shown  in  Fig. 2a.  It  indicates  a property  of  F(X) 
which  follows  from  our  requirements:  For  |x|  beyond  some  selected  non-negative 
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We  subsequently  consider  cases  of  functions  which  are  differentiable  but 
not  analytic  at  X = 0. 
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value  which  we  will  call  X (no  relation  to  the  Lame  constant)  which  is  shown 

o 

in  Fig. 2a,  we  may  regard  |F(X)|as  approximately  constant  and  equal  to  1.  We  may 
select  as  large  as  we  deem  necessary  in  order  to  consider  F(X)  constant 
for  |x|  > X . 

We  define  a quality  factor  Q by  the  relation  Q = 1 - F(X^).  Then  Q is  the 
upper  bound  of  relative  error  from  F being  exactly  equal  to  one  which  we  are 
willing  to  accept.  Exactly  constant  behavior  of  F(X),  corresponding  to  perfect 
quality,  would  have  Q = 0;  while  a total  lack  of  quality,  which  would  regard  any 
deviation  of  F(X)  from  unity  as  negligible  and  which  would  enable  us  to  regard 
any  F(X)  value  on  [0,1]  as  constant  in  X and  equal  to  1,  would  correspond  to 
Q = 1.0  or  100%.  Within  the  scope  of  our  theory  we  may  choose  any  Q within  the 

range  0 < Q < 1,  with  Q = 0 not  an  admissible  choice.  Thus  we  are  free  to  select 

the  quality  of  our  approximation  to  constant  behavior  (constant  slope)  over  a 
continuous  range  from  total  lack  of  quality  to  any  arbitrary  nearly  perfect  but 
never  exactly  perfect  quality. 

We  may  select  Q = .5%  if  we  wish  to  regard  F(X)  as  constant  and  equal  to 
one  when  its  true  magnitude  is  on  the  interval  (.995,  1.0]  and  F(X)  as  noncon- 
stant, varying  with  X,  when  its  magnitude  is  less  than  .995;  this  clioice  of 
Q = .5%  would  then  correspond  to  a choice  of  X^  = P ^ (1  - Q)  = F ^(.995). 

We  see  that  the  choice  of  Q,  with  Q restricted  to  0 < Q ^ 1,  and  the 

particular  function  chosen  for  F determine  uniquely  the  value  of  X^.  Henceforth 
we  will  use  the  term  "constant  behavior"  to  mean  behavior  which  is  approximately 
constant  to  an  accuracy  commensurate  with  the  value  of  the  factor  Q,  and  we  will 
use  the  term  "linear  behavior  in  X"  to  denote  the  behavior  of  a function  which  is 
approximately  linear  in  X,  and  where  the  slope  of  this  function  is  tlnen  approxi- 
mately constant  to  a degree  of  accuracy  coTumensurate  with  the  value  of  the 
quality  factor  Q. 

Appendix  A lists  particular  F(X)  functions  which  we  have  used  so  far. 
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Kernel  Forms 

We  now  introduce  two  sets  of  kernel  functions  of  X:  the  set  U^(X),  V^(X) 
and  the  set  U2(X),  ^2^^^  which  enable  us  to  develop  two  distinct  forms  of  Y(X). 

The  First  Set  of  Kernel  Forms 


Define 


X + X 

"i »'  - K 

f 1 


X.  - X 


'i »'  - *o  + "‘o  (xT^) 


Note  that  U^(0)  = V^(0)  > 0 since  ^ ^ 0-  Recalling  the  form  of  F(X),  we 

have  for  F(Uj^(X))  and  F(Vj^(X))  the  forms  shown  respectively  in  Figs.  3a  and  4b. 
We  may  interpret  U^(X)  as  a translation  of  the  vertical  axis  and  a rescaling 
(stretch)  of  the  horizontal  (X)-axis;  and  we  may  interpret  as  a translation 
of  the  vertical  axis  and  a reflection  and  rescaling  (stretch)  of  the  horizontal 
(X)-axis.  We  now  add  the  two  functions  F(U^^(X))  plus  F(V^(X))  to  obtain  the 
form  indicated  in  Fig. 3c.  From  the  required  behavior  of  the  F-functions, 
and  as  is  indicated  in  Fig.  3c,  the  sum  (F(U^)  + F(V^))  is  constant  on  tlie 
interval  |x|  < X^  and  again  constant  on  the  unlaounded  interval  |x|  > X^  to 
within  an  accuracy  commensurate  with  the  choice  of  quality  factor  Q.  It  is 
extremely  important  to  observe  that  while  these  two  intervals  of  constant 
behavior  of  F(X)  must  be  joined  by  an  interval  of  nonlinear  behavior  whose 
curvature  depends  upon  the  particular  choice  of  the  function  F,  the  intervals 
of  constant  behavior  may  extend  beyond  the  value  of  ±X^  and  may  commence 
before  |x|  attains  the  value  X^.  That  is,  X^  and  X^  absolutely  do  not  denote 
the  points  at  which  F(X)  respectively  begins  and  ends  nonconstant,  nonlinear 
behavior;  rather  they  specifically  denote  those  |x|  points  at  which  F(X)  must 
for  |x|  values  less  than  X^  or  larger  than  X^  be  constant  in  behavior.  The 
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particular  choice  of  F will  determine  the  size  of  any  excess  interval  of  constant 
behavior  in  |x|  exceeding  or  preceding  X^. 

We  now  obtain  V (X)  by  multiplying  (F(U^)  + F(V^))  by  a scale  constant 
and  adding  the  constant  to  this  sum: 

Y*  (X)  = + A^(F(U^(X))  + F(V^(X))  . (3) 


Since  (0)  = (0)  and  since  Y'  (0)  = we  have 

(S,  - S^) 

\ 2F(Uj^(0))  • 


Clearly  Y'  (X)  is  constant  over  the  desired  intervals  since  Y'  (X)  is  con- 
structed from  the  sum  (F(U^)  + P(V^))  which  is  constant  on  these  intervals. 
Let  T(X)  denote  the  integral  of  F(X) 


T(X)  = F(5)  d5  . 


Then  from  (3)  we  have 


Y(X)  = S^X  + B^(T(U^(X))  - T(V^(X))  ) 

Bl  = - X.)/(2X^) 


with  0 < Q 5 1,  = f“  (1-Q),  and  with  ^ > 0.  as  defined  in 

Fig.l. 

Although  the  parameters  X^,  X^,  S^,  may  be  selected  to  be  material 
parameters  which  characterize  a given  specific  material,  X^  is  absolutely  not 
a material  parameter.  The  parameter  X^  scales  the  kernels  of  the  Y' -function 
to  obtain  the  desired  constant  behavior,  and  consequently  X^  is  a quality- 
controlling variable  which  must  be  selected.  A change  in  the  selection  of 
the  F-function  used  to  represent  Y' (X)  and  Y(X),  for  a fixed  chosen  quality 


; 1 » i ' 


factor  will  result  in  a change  in  the  value  of  X^. 


For  the  special  case  of  some  limiting  values  of  the  parameters  X^,  S^, 
S^,  Y(X)  in  (5)  can  yield  indeterminate  limits.  In  this  case  we  reconstruct 
Y(X)  from  the  limiting  parameter  form  of  Y' (X)  in  (3)  which  is  then  integratea. 
The  results  are  given  in  Appendix  B. 

The  Second  Set  of  Kernel  Forms 

An  alternative,  useful  set  of  kernel  forms  is 

U2(X)  =-  \^  + R(X^  + X)  (7) 

and 

V- (X)  = - X + R(X  - X)  (8) 

2 of 


where  R > \^/X^.  Here  we  have  completely  abandoned  utilization  of  the 
parameter  X^  and  have  instead  elected  to  introduce  a new  parameter  R,  not 
related  to  X^,  which  we  shall  refer  to  as  the  "amplitude  constant"  and  which 
is  subject  to  the  constraint  R > X^X^^.  Figures  2b  and  2c  illustrate  graphi- 
cally the  forms  of  F(U2(X))  and  F(V2(X)).  As  in  the  case  of  the  previous 
kernel  forms,  we  may  regard  the  kernel  as  a translation  of  the  vertical  axis 
and  a rescaling  (stretch)  of  the  horizontal  (X)-axis,  and  we  may  regard  as 
a translation  of  the  vertical  axis  and  a reflection  and  rescaling  (stretch)  of 
the  horizontal  axis.  Again  we  add  two  F-function  forms  to  obtain  the  sum 
(Fdl^)  + FCV^))  which  appears  graphically  in  Fig.  2d.  Note  that  in  this  case,  as 
a result  of  the  use  of  our  new  kernel  choices,  although  the  resultant  function 
sum  is  again  initially  constant  on  a symmetric  interval  about  the  origin,  we  can  I 

no  longer  discuss  or  suggest  a minimum  bound  (X.)  upon  the  size  of  this  interval.  ^ 
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However,  since  we  have  retained  use  of  the  parameter  in  our  new  kernels,  and 
since  R is  greater  than  since  F(X)  is  constant  for  |x|  greater 

than  we  are  able  to  conclude,  exactly  as  in  the  case  of  our  previous  kernel 
forms,  that  the  function  sum  (F(U2)  + F(V2))  is  constant  for  all  X on  tlie 
unbounded  intervals  |x|  > X^,  Multiplying  the  sum  (F(U2)  + FCV^))  by  a scale 
constant  and  adding  the  constant  enables  us  to  represent  Y' (X)  in  terms 
of  the  second  set  of  kernels: 


y (X)  = S,  + (F (U. (X) ) + F (V- (X) ) ) 


since  Y'  (0)  = and  > 0 we  have 


S.  - S- 

A = ^ ^ 

2 2F(U2(0)) 


and  integration  of  (9)  in  terms  of  T(X)  results  in 


Y(X)  = S^X  + B2(T(U2(X))  - TCV^CX))  ) 


= A^/R  . (12) 

Clearly  Y(X)  is  linear  in  X for  |x|  > X^  because  of  our  synthesis  of  Y(X) 
from  the  function  sum  (F(V^)  + F (V^) ) which  is  constant  on  the  interval 
|x|  > X^;  because  of  this  linearity  we  will  refer  to  the  region  |x[  > as 
the  region  of  steady-state  behavior  of  Y{X),  or  the  steady-state  region. 
Because  Y(X)  is  linear  in  its  steady-state  behavior,  the  effect  of  varying  the 
selection  of  the  amplitude  parameter  R is  to  translate  vertically  the  steady- 
state  Y-curve.  Then  over  the  range  of  admissible  R and  the  limiting  R-va)ues. 
i.e.,  Rs  tX^/X^,«>]  , for  any  fixed  choice  of  S^,  S^,  positive  X^,  and  steady- 
state  X,  there  is  one  R-value,  denoted  as  R . , which  renders  the  smallest 
possible  value  of  |y(X)|,  and  there  is  one  R-value.  denoted  as  R . which 
renders  the  largest  possible  value  of  steady-state  amplitude  |y(X)|.  The 
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values  of  R . and  R are  independent  of  the  specific  values  of  S.  and  S- 
min  max  if 

and  depend  only  upon  the  ranges  of  and  S^.  Examination  of  the  Y (X ) -function 

indicates  that  for  the  cases  positive  and  greater  than  or  negative 

and  less  than  S^,  R . is  finite  and  R = a>;  for  the  converse  cases  of  S. 

f’  min  max  i 

positive  but  less  than  or  negative  and  greater  than  S^,  is  finite 

and  R . = oo.  Consequently,  the  effect  of  the  choice  of  S . , S-  combinations 
min  ^ , i’  f 

determines  whether  R = ;»  corresponds  to  R or  to  R . ; and  whether  a particular 

^ max  min 

finite  critical  R-value  corresponds  to  R , or  to  R 

min  max 

Throughout  this  report  we  will  consider  the  case  5 0.  In  this 

case  R =00  and  R . must  either  be  the  lower  limit  of  the  range  of  admis- 
max  min 

sible  R,  equal  to  X /X,.,  or  R . must  be  an  R-value  which  satisfies  ' .le  on- 
’ of’  min 

straint  dY(X;R)/dR  = 0. 

Since  the  numerical  value  of  is  independent  of  the  particular  values 

of  and  while  depending  upon  the  ranges  of  these  values,  the  numerical 

value  of  then  depends  only  upon  the  F-function  used,  the  chosen  quality 

factor,  the  value  of  X^,  and  the  choice  of  the  steady-state  X-value.  As  a 

consequence  of  the  linear  behavior  of  Y(X)  in  its  steady-state  regions,  varying 

the  choice  of  X within  the  steady-state  region  will  not  change  the  minimizing 

value  of  R . . By  considering  the  limiting  case  of  X ->  co  we  obtain  a trans- 
min 

cendental  relationship  enabling  us  to  write  and  use  R . = R . (X_)  for  any 

min  min  f 

given  particular  F-function  choice  and  desired  quality  factor  Q.  A derivation 

of  the  relation  R . (X,)  for  a particular  function  choice  appears  in 

min  f 

Appendix  C. 

While  a decrease  of  R below  the  value  R . and  on  the  interval  R . > 

min  min 

R S result  in  an  increase  in  the  steady-state  amplitude,  this  range 

of  amplitude  is  contained  as  a subset  of  the  amplitude  range  corresponding  to  R 
increasing  beyond  Consequently  for  X in  the  steady-state  regions  and 
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fixed,  increasing  R beyond  allows  us  to  sweep  the  entire  range  of 

available  steady-state  amplitudes.  The  parametric  dependence  of  Y(X)  upon  the 

amplitude  parameter  R is  illustrated  in  Figs. 4 and  5.  In  Fig. 4 the  band  of 

available  Y(X)-curves  for  R . s R oo  is  shown.  Figure  5 shows  the  variation 

min 

of  Y as  a function  of  R and  the  value  of  R , is  clearly  identifiable.  It  is 

min 

further  evident  that  Re [R  . , contains  the  amplitude  values  of  steady- 

min 

stat  ' Y(X)  for  Re  [1  /X.,  R . ]. 

o f’  min 

Because  of  this  particular  parametric  dependence  of  the  Y (X) -function 

upon  R,  we  now  propose  to  reduce  our  admissible  range  of  R to  Re [R  . , ; no 

min 

other  R- values  need  be  considered. 
biscussion 

Both  sets  of  kernel  forms  can  be  used  with  any  of  the  functions  listed  in 
Appendix  A.  A great  variety  of  F(X)  function  choices  are  at  our  disposal.  If 
necessary,  the  set  of  available  function  choices  can  be  enlarged  by  adding  new 
functions  satisfying  the  requirements  placed  upon  the  form  of  F(X),  and  by 
considering  linear  combinations  of  available  F(X)  forms. 

To  describe  a particular  Y(X)  four  independent  parameters  must  be  selected. 
The  initial  and  final  slope,  and  S^,  respectively  are  common  to  both  kernel 
sets.  In  addition,  for  the  first  kernel  set  X^  and  X^  can  be  prescribed  pro- 
vided X-  S X.  >0.  In  the  case  of  the  second  kernel  form  R and  X^  can  bo 
f 1 f 

selected  subject  to  the  restriction  Re [R  . , “]  where  R . S \ /X,. 

min  min  o f 

Both  kernel  forms  yield  linear  (B2, Bll, B20, B24 ) or  bilinear  behavior  (B15 
and  B28)  for  limiting  values  of  the  parameters.  These  "linear"  forms  of  Y(X) 
are  independent  of  F(X). 

To  obtain  numerical  values  of  Y(X)  calculators  must  be  used.  For  the 
hyperbolic  tangent,  the  inverse  tangent  and  the  rational  function  representations 
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It 
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of  F(X),  Y' (X)  and  Y(X)  can  be  easily  calculated  using  an  "electronic  slide-  j 

rule".  We  have  found  it  useful  to  employ  a central  computer  and  illustrate  j 

the  results  by  graphs  of  Y' (X)  and  Y(X).  Some  of  them  are  reproduced  below  to  | 

illustrate  the  properties  of  the  functions. 

'j 

Figure  6 shows  the  behavior  of  the  normalized  Y'  (X)  function.  Eg.  (3),  for  | 

F(X)  based  on  the  error  function  (Al)  and  the  hyperbolic  tangent  function  (AS), 

respectively.  The  constant  behavior  for  |x|/X^  < 1 and  |x|/X^  > 3 for  both 
functions  is  evident  as  well  as  the  different  behavior  in  the  nonconstant  i 

region.  (The  ratio  p,  = X^/X^  determines  the  interval  of  constant  behavior.  | 

We  emphasize  that  Q = 0.005  in  both  cases.) 

The  influence  of  the  selection  of  u,  = X../X,  on  the  behavior  of  Y'  (X)  is  [ 

illustrated  in  Fig.  7.  As  |j,  increases  from  1.1  to  3.5  the  transition  from  S.  t 

to  becomes  more  gradual.  This  graph  also  shows  that  although  constant  I 

behavior  is  definitely  obtained  for  [x(/X^  < 1 and  |x(/X^  > it  may  extend  | 

beyond  jx|/X^  = 1 and  commence  before  |x|/X^  = p,.  (The  comparison  between  j 

the  p.  = 1.1  and  the  p = 3.5  curves  illustrates  this  point  very  well.)  Although  j 

'1 

the  transition  for  p = 1.1  appears  very  sharp  on  the  graph,  the  Y' (X)  for  J 

p = 1.1  is  analytic.  (Note,  that  a jump  discontinuity  would  be  obtained  at  1 

■f 

|x|/x^  = 1 for  p = 1,  Eg.  (B14).  S 

Figure  8 shows  the  parametric  dependence  of  Y' (X)  on  the  ratio  S /S.  | 

illustrating  that  and  may  be  arbitrarily  selected.  i| 

The  next  three  graphs  are  illustrations  of  tlie  capabilities  of  the  second  | 

set  of  )cernel  forms  for  Y(X).  In  these  illustrations  Y(X)  is  identified  with 

Cauchy  stress  ct  and  X assumes  the  role  of  strain  e because  of  the  obvious 

resemblance  of  the  Y(X)  vs  X-curves  to  the  stress-strain  diagram.  j 

Figure  9 shows  a/S . vs  strain  based  on  Egs.  (11)  and  (Al)  for  various  values  ij 

^ [! 

of  X„  for  fixed  S.  and  S^.  At  each  value  of  X.,  R = R . , so  that  we  have  J 

f 1 f f’  min’  j| 
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drawn  the  lowest  possible  curve  for  each  X^-value.  Therefore,  at  each  value 

of  higher  curves  could  be  obtained  by  increasing  R beyond  see  Fig. 4. 

For  X^  = 0.002  we  have  graphed  the  bilinear  limiting  form  for  R = oo  in 

addition  to  the  curve  for  By  increasing  R beyond  Y(X)-curves 

higher  than  the  curve  for  R . but  lower  than  the  curve  for  R = can  be 

min 

obtained. 

In  the  discussion  of  the  rational  function  form  of  Y(X),  Eg.  (A16 ) , we 

mentioned  that  each  choice  of  the  exponent  5 determines  a different  function. 

This  is  illustrated  in  Fig. 10  where  stress-strain  curves  for  f ==  3 and  f = 11 

for  the  same  X_-values  are  displayed,  respectively.  At  each  X^,  R=R  . (X^). 

t f’  min  r 

It  is  seen  that  the  curves  for  5 = 11  are  higher  than  those  for  5=3.  By 
coincidence  the  curve  for  X^  = 0.004  and  5 ~ ^ overlaps  almost  everywhere  with 
the  curve  X^  = 0.001  and  5 = H-  The  curve  for  X^  = 0.001  and  R = 00  is  drawn 
for  comparison.  This  limiting  curve  applies  to  both  values  of  the  exponent  f. 
The  flexibility  of  the  rational  function  is  obvious. 

Finally  Fig. 11  uses  a different  scale  for  the  strain  and  a negative  value 
for  and  illustrates  that  this  choice  is  also  possible  for  Y(X)  based  on 
Eg.  (11). 

Application 

From  the  discussion  of  Figs. 9 -11  it  is  obvious  that  the  Y (X) -functions 
developed  from  the  second  set  of  kernel  forms  are  most  suitable  for  approxi- 
mating experimentally  determined  stress-strain  diagrams.  The  slopes  S^, 
as  well  as  the  value  of  X^  are  immediately  obtained  from  these  diagrams.  For 
each  F (X)-function  can  be  calculated  by  the  procedure  given  in  Appendix  C. 

The  resulting  minimum  Y(X)-curve  can  be  computed  and  compared  with  the  diagram 

to  be  matched.  If  the  diagram  is  above  Y(X)  for  R = R . , the  R that  results 

min’ 


in  the  best  approximation  can  be  found.  If  the  diagram  is  below  Y(X)  for 

R = R . , another  function  F(X)  must  be  chosen  and  the  procedure  repeated. 

min 

We  recall  that  the  use  of  the  second  set  of  kernel  forms  does  not  permit 
an  easy  identification  of  the  extent  of  the  linear  behavior  of  Y(X)  in  the 
vicinity  of  X = 0.  It  would  therefore  appear  that  the  first  kernel  form  would 
offer  advantages  if  an  initial  interval  of  almost  linear  behavior  is  of  interest. 
However,  the  choice  of  the  first  kernel  set  permits  the  specification  of  this 
interval  but  does  not  allow  us  to  control  the  amplitude  of  the  steady-state 
Y(X)-curve.  It  is  for  this  reason  that  the  second  set  of  kernel  forms  is 
preferred  for  application  to  stress-strain  diagrams. 

Another  application  of  the  Y' (X)  and  Y (X) -functions  is  in  the  identifica- 
tion of  coefficient  functions  of  nonlinear  constitutive  laws  such  as  the  dif- 
ferential form  proposed  in  [4],  Eg.  (3).  Ihe  Y (X)-function  is  then  identified 
with  g (e ) and  the  Y' (X) -function  assumes  the  role  of  m(a).  This  function  is 
responsible  for  modeling  the  nonlinear  stress  dependence  of  the  creep  rate 
which  is  observed  in  metals,  see  Eq.  (9)  in  [4] . It  is  shown  in  [4]  that 
"numerical  experiments"  using  the  same  functions  but  different  values  of 
parameters  can  change  the  "character"  of  the  results.  As  a specific  example, 

= 0 in  the  Y(X)  used  for  g (e ) reproduces  secondary  creep  whereas  > 0 
produces  primary  creep  only;  see  Eq.  (10)  and  Figs. 7 and  8 in  [4]. 

Also  the  almost  linear-elastic,  rate-independent  initial  range  found  in 
stress-strain  diagrams  obtained  at  various  strain  rates  can  be  reproduced  with 
the  differential  form,  Eq.  (3)  in  [4]  if  g (e ) has  the  initial  almost  linear 
region  of  Y(X)  and  m (ct)  and  k(e)  are  almost  constant  as  is  the  case  for  the 
Y' (X) -functions ; see  Figs.l  and  2 in  [4]. 

A further  potential  application  ’s  in  the  modeling  of  history  dependence 
in  the  sense  of  plasticity.  By  making  X^  and  R depend  upon  prior  deformation 


history,  which  is  possible  through  the  structure  memory  function  the 


graphs  in  Fig. 9 could  be  interpreted  as  stress-strain  diagrams  after  various 


amounts  of  prior  cold  work 


ture  and  the  graphs  in  Fig. 9 can  be  interpreted  as  having  been  obtained  at 


various  temperatures.  Although  S.  and  S are  kept  constant  in  Fig. 9,  they  can 


be  changed  to  any  desired  value  if  such  a need  arises  in  application 


(viscoplasticity) . Their  closed  form  makes  them  comparatively  easy  to  use 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Figure  7 


Figure  8 


Figure  9 


Schematic  representation  of  the  function  Y(X)  and  iti^;  derivative 
y (X)  (stress-strain  diagram).  For  |x(  < X^  and  for  |x|  > X^ 

almost  linear  (constant)  behavior  is  observed. 

Schematic  illustration  of  F(X),  Fig. 2a;  the  axes  are  translated 
stretched  (Fig. 2b);  translated,  reflected  and  stretched  (Fig. 2c). 

The  addition  of  F (U^ ) and  F (V^ ) yields  within  a constant  the 

desired  Y'(X)-form,  Pig. 2d. 

The  axes  are  translated  and  stretched  by  (3a)  and 
translated,  reflected  and  stretched  by  (3b)  before  F(U^) 
and  F(V^)  are  added  (3c).  Figure  3c  represents  (within  a 
constant)  the  desired  Y' (X)-form. 

Schematic  showing  the  dependence  of  Y(X)  on  the  parameter  R for 

R . ^ R s CO.  Note  that  for  R = ” bilinear  behavior  is  obtained, 

min 

The  other  parameters  S^,  S^  and  X^  are  held  constant.  ; 

The  dependence  of  Y on  the  parameter  R for  X = .003  at  three  j 

different  values  of  X^.  The  curves  are  based  on  F(X)  = ERF(X)  ! 

and  exhibit  the  minimizing  value  of  R.  In  our  application  we 

restrict  R to  R > X /X,.  1 

° ^ 

The  normalized  function  Y'  (X),  Eq.  (3),  based  on  the  error  function  J 

Y^^^,  Eq.  (A3),  and  on  the  hyperbolic  tangent  function  1 

Eq.  (A7);  M'  = X^/X^  = 3,  S^/S^  = .25,  Q = 0.005  in  both  cases.  i 

The  influence  of  the  choice  of  ti  = X^/X^  on  Y'  (X).  Y'  (X)  is 
based  on  Eqs.  (3),  (A3)  and  S^/S^  = .25,  Q = 0.005. 

Y'  (X)  at  different  values  of  S -/S  . with  con.stant  p = X-/X.  =3.  ! 

f 1 f 1 

Y'  (X)  is  based  on  the  error  function,  Eqs.  (A3)  and  (3),  Q = 0.005. 

o/S. -curves  based  on  Eqs.  (11)  and  (A4).  At  each  X^.  R = R , (X,.).  I 

1 f min  f ' 

For  X = 0.002  the  bilinear  form  (R  = ®)  is  also  shown, 
a ^ Y(X);  X ^ e (strain)  ; S^  = 21.7  X 10  psi  (=*  150  GPa), 

S^/S.  = 0.02,  Q = 0.005. 

f 1 • . 


Figure  10 


Figure  11 


a/s ^-curves  based  on  Eqs.  (11)  and  (A16 ) for  § = 3 and  F = 11  at 

various  X_.  For  each  X.,  R = R . (X  ) . For  X,  = .001  the 
f f’  min  f f ^ 

bilinear  form  (R  = “>)  is  also  shown.  = 21.7  X 10  psi 

(<=«  150  GPa),  S^/S^  = 0.02,  Q = 0.005. 

a/S^-curves  based  on  Eqs. (11)  and  (A4)  for  negative  S-values 

and  expanded  scale.  R = R . (X-).  For  X^  = 0.01  the  bilinear 

form  (R  = o>)  is  also  drawn  in.  = 21.7  X 10  psi  150  GPa), 

S^/S.  = -0.005,  Q = 0.005. 
f 1 ’ 


iiraa* 


APPENDIX  A 


PARTICULAR  CHOICES  FOR  THE  F(X),  T(X),  Y' (X)  AND  Y(X)  FUNCTIONS 

For  convenience  wp  will  be  using  generalized  kernels  U,  V and  constants 


A,  B.  Iben  if  use  of  the  first  kernel  set  is  desired,  the  set  U,  V,  A,  B are 
to  assume  the  values  of  the  set  U^^,  A^,  B^  as  defined  in  Eqs.  (1),  (2),  (4), 

(6).  For  the  second  kernel  set  U,  V,  A,  B will  assume  the  values  given  in 
Eqs. (7),  (8),  (10),  (12). 


1)  Error  Function 


F(X)  = erf(X)  with  erf  (oo)  = l 

(Al) 

1 2 

T(X)=Xerf(X)  + exp(-  X ) 

yir 

(A2)  ; 

Y*  (X)  = + A(erf (U)  + erf  (V)) 

(A3)  \ 

; 

1 

Y(X)  = S.X  + B(U  erf(U)  - V erf  (V)  +— ^ (exp(-U^)  -exp(-V^))) 

i 

(A4) 

2)  Hyperbolic  Tangent  Function 

F(x)  = tanh(X) 

(A5) 

T(X)  = S/n  (cosh(X))  is  base  e) 

(A6) 

Y'  (X)  = S^  + A (tanh (U)  + tanh  (V) ) 

(A7) 

Y(X)  - S,X  + B 

f \cosh(V)/ 

(A8) 

4)  Rational  Function 

F(X)  = sign(X)(^l  - (l  + |x|)  ) , I > 2 (A! 

T(x)  = |x|  + 7|Tiy  (1  + (a: 

Y'  (X)  = + A [sign(U)  (1  - (1  + |ul  )^”^)  + sign(V;  (1  - (1  + | v|  "’)] 


Y(X)  =S-X+  ■■  [(5-2)(|u|  - lv|)  + (1  + |U|)^~^  - (l+lv|)^"^] 

f (^-2) 


5)  Incomplete  Gamma  Function 

F(X)  = sign(X)  Y , |xr)/r(i),  0)  > 0 


1 [ 

1 1 A 

1 1 / 2 

1 1 ud\"| 

X yI“  , 

' \au  ’ 

i>'i  j - ns  > 

i'=i  ;J 

Y'  (X)  = + -7^  [sign(U)  , | U | ‘“^  + sign  (V)  Y (^  , | V | (A19) 


Y(X)  = S^X 


_ r 

1 1 /I 

1 

u y(-  , 

u j 

' \UD  ’ 

- iviv(^  ivr; 


(I  - 1-1“^)  ^-(1^ 


where  Y is  the  incomplete  gamma  function  [5]  and  sign(X)  is  defined  as 

1 X > 0 


sign(X)  = { 0 X = 0 


-1  X < 0 


Here  we  introduce  F (X) -functions  which  are  not  analytic  but  are  differentiable 
at  X = 0. 


I • i 


For  both  the  rational  and  incomplete  gamma  function  each  value  of  § 
and  u),  respectively,  represent  one  possible  y(X)-form.  Therefore  an  infinite 
class  of  functions  with  different  nonlinear  behavior  is  available. 


'I  ?*►.. 


APPENDIX  B 

y (X)  and  Y(X)  FUNCTION  REPRESENTATIONS  FOR  LIMITING 
VALUES  OF  THE  PARAMETERS 


1 


In  some  cases  the  limiting  parameter  values  result  in  indeterminate 
forms  of  Y(X)  when  Eqs. (5)  or  (11)  are  used  to  represent  Y(X).  This  is  a 
consequence  of  the  construction  of  Y(X)  in  (5)  or  (11)  from  Y' (X)  in  (3) 
or  (9)  for  nonlimiting  pcurcimeter  values.  Consequently,  when  these  indetermi- 
nate limits  arise,  we  reformulate  Y(X)  as  the  integral  of  the  limiting 
parameter  form  of  Y' (X)  in  (3)  or  (9);  this  procedure  is  fundamentally  con- 
sistent with  our  construction  of  Y(X)  from  Y' (X).  Alternatively  the  inde- 
terminate limiting  form  of  Y(X)  may  be  resolved  when  specific  forms  of  Y(X) 
are  chosen;  and  in  some,  but  not  all  cases,  general  forms  of  Y(X)  can  be 
used  for  the  resolution  of  the  indeterminate  limits. 


Both  sets  of  Icernel  forms. 


Limit  -*  S . . 

f 1 


From  Eqs. (3),  (4),  (5),  (6)  and  (9),  (10),  (11),  (12),  we  obtain 


lim  Y' (X)  = S . 

S^-*S.  ^ 

f 1 


(Bl) 


lim  Y(X)  = S.X 

Sf-Si 


(B2) 


First  set  of  )cernel  forms. 


Limit  X . -*0. 
1 


From  (1),  (2),  (4),  (6),  the  following  limits  are  easily  derived 


lim  U (X) 

X,  -0  ^ 

1 


X + 2\  X/X^ 
o of 


(B3) 


lim  V (X) 

X.-O  ^ 

1 


X ~ 2\  X/X, 
o of 


(B4) 


\-i 


t 


Then  Y 


where  the  above  limiting  kernel  functions  and  parameter  values  are  used 


Limit  X 


The  following  limits  result 


Proceeding  consistent  with  our  construction  of  Y(X)  as  the  integral  of  Y' (X) 


we  integrate  (BlO)  and  obtain 


Limit  X 


for  X > X.  lim  U (X)=+®;  lim  V (X)  = 

^ X.-X^  ^ X.-X.  ^ 

if  if 


for  |x|  < X^  lim  U,  (X)=+a>;  lim  V (X)=+< 

' ' f X.-X^  1 X.-X.  1 

if  if 


for  X < -X,  lim  U.(X)=-®;  lim  V (X)=+' 
^ X.-X^  X.-X 

if  If 


lim  A = (S.  - S )/2 
X.-X.^  lit 
1 f 


Proceeding  as  before,  we  obtain  first 


lim  Y'  (X) 

X.-X. 

1 f 


S.  X < X, 


■Sf  |X|>X^ 


and  integrate  to  obtain 


lim  Y(X) 

X.-X^ 

1 f 


S.X^  + S^(X-  X.) 
if  f f 


-S.Xf  + S^(X  + X^) 


* f 

|x|  ^ x^ 

X S -X, 


It  is  significant  to  note  that  bilinear  behavior  is  obtained  in  the  limit 
X^  - X^  independent  of  both  the  particular  F (X) -function  and  the  quality 
factor  Q. 


Second  set  of  kernel  forms 


Limit  X^  — 00. 


FromEqs.(7),  (8),  (9),  (10),  (11)  we  obtain 


lim  U (X)  = + 00 

X^-oo  2 


lim  V (X)  = + 00 
X^-oo  2 


(B17) 


B4. 


k 


lim  A_  = (S.  - S,)/2 

X -CO  2 '1  f 

f 

(B18) 

lim  Y' (X)  = S. 

X^-00  1 

(B19) 

lim  Y(X)  = S.X  . 

(B20) 

i 1 

X^-00  1 

[ 

f 

l 

I 

1 

i 

1 


I 

I- 


We  note  then  that  for  each  set  of  kernels,  as  “ Y(X)  approaches  a func- 

tion exactly  linear  for  all  X,  with  slope  S^. 

Limit  X^  -•  0. 

For  the  second  kernel  set,  when  X^  -•  0,  we  must  recall  the  restric- 
tion R > \ /X  ; that  is,  while  X,  -•  0 the  product  RX.  attains  a finite  value 
o f ’ f t 

greater  than  the  constant  Consequently  R must  approach  infinity  as  X^ 

approaches  zero.  Then  we  have 


lim 

U„  (X)  = +00  ; lim  V, 

X^-*0 

2 0 

R-.00 

R-*oo 

1 

lim 

U~  (X)  = - 00  ; lim  V, 

Xf-0 

^ X^-  0 

R-00 

R-.00 

lim 

0 ^ 

. . * 
= a positive  constant 

R-.00 

(B21) 


(B22) 


Because  of  the  lindts  involved,  the  particular  value  of  this  constant  will 
not  enter  into  Y'  or  Y.  Further,  note  that  is  fundamentally  defined  using 

U„(0)  and  consequently  in  performing  the  limit  of  A we  must  use  lim  U (0) 

2 2 0 

R-.0O 

rather  than  using  limits  on  U2  then  letting  X — 0. 


Consequently 


and 


lim  y (X) 


R-.00 


S 


f 


Lim  Y(X)  = 
X^-  0 

pt-oo 


SfX 


Limit  R -*  00,  X,  finite,  positive 


We  have 

for 

X 

^^f 

lim 

R-OO 

for 

|X| 

I 

lim 

R-»“ 

for 

X 

<-^f 

lim 

R-*oo 

and 

lim 

R-.00 


u.  (X)  = +“  ; 

lim 

V,  (X) 

= ■»  00 

2 

R-oo 

2 

U.  (X)  = +“  ; 

lim 

V_  (X) 

= + 03 

2 

R-.00 

2 

U.  (X)  =-«  ; 

lim 

V-  (X) 

= +C0 

2 

R-*00 

2 

^2  " ^^i  ■ • 


Therefore 


and 


r ^ ■ 

lim  y (X)  = j ^ 
R-»oo  ^ S 


lx|  < x^ 
|x|  > x^ 


lim  Y(X) 

R-.0O 


S.X-  + S,(X  - X.) 
1 f f f 

i ^i^ 

-S.X,  + S-(X  + X.) 
iff  f 


X ^ 

lx|  ^ x^ 

X i -x^  . 


(B23) 


(B24) 


(B25) 


(B26) 


(B27) 


(B28) 


We  note  then  that  taking  the  limit  as  R approaches  infinity,  results  in  a 
bilinear  form  for  the  Y(X) -function  and  Y'  (X)  is  discontinuous  at  X = X^. 
The  maximum  amplitude  |y(X)|  at  any  steady  state  range  of  |x|,  for  fixed  X 


B6. 


is  then  S^|x|  + (S^  - S^)  corresponding  to  the  steady-state  amplitude  of 

...  * 
the  bilinear  R = ® response  . 

This  bilinear  Y(X)  response  for  the  U^,  kernel  representation  and  the 
limit  R = 00  is  identical  to  the  bilinear  Y(X)  response  which  results  when  the 
^1’  ^1  used  with  the  limit  X^  • 

Limit  R - X /X,,  X_  finite,  positive 
o f’  r 


lim  U_(X)  = X X/X, 
^ ,..2  of 


lim  V_ (X)  = - X X/X, 
R-X  /X^  2 of 

o f 


(B29) 


(B30) 


lim  A = 00  (b31) 

R-*X  /X,  2 
o f 


lim  (F(U.)  + F(V-))  =0  (B32) 

R^X  /X^  2 2 

o f 


In  this  limit  both  our  Y'  and  Y-forms  in  (9)  and  (11)  are  indeterminate. 


Resolving  this  limit  we  have 

F'  (X  X/X  ) 

lim  YMX)  - S,  ♦ (S,  - S,)  ..  (B33) 

n-yXf 


lim  Y(X) 


K'J) 


S^X  + 


(Si  - s^) 


X,,F(X  X/X,) 
f o f 


X F' 
o 


(0) 


(B34) 


t 

f 


For  S.  > S.  a 0. 
1 f 


APPENDIX  C 


DEVELOPMENT  OF  THE  REPRESENTATION  R . (X^) 

mm  f 


Herein  we  will  construct  the  representation  of  the  X -dependence  of 


for  the  case  v/hen  R . is  finite  and  R 
mm  ; 


max 


mm 


complementai-y  case  of  either  2 0 cr  s 0.  the  values  of  R 


min 


is  the  finite  R-critical  point 


mm 


max 


max 


Note  that  R . denotes  the  R-value  that  minimizes  the  steady-state  amplitude 
min 

of  Y(X)  and  denotes  the  corresponding  maximizing  value  of  R. 


Proceeding,  we  regard  Y as  a function  of  R,  write  Y(X;R)  and  restrict 


X to  |xj  > X ; then  for  R > \ /X_,  Y is  an  analytic  function  of  R 


we  seek  a finite  R > \ /X-  which  is  a solution  to 


min 


and  which  satisfies 


From  Eq.  (Cl)  and  Eqs.  (7)  - (12)  we  have,  for  R > X /X 


F(U-)(X-+X)  - F(V_)(X.-X)  - (1/R)(T(U.)  - T(V_))(1  + 


is  independent  of  our  particular  choice  of  steady-state  X 


min 


Then  to  greatly  simplify  the  relation  (C3)  and  the  representation  of  (C2) 


This  is  the  case  for  all  functions  proposed  in  Appendix  A 


C2. 


we  perform  the  limit  of  these  relations  for  X -*  oo  and  introduce  the  para- 
meters B and  B defined  as 


3 = RX^  - X 20 
f o 


3 = R . X^  - X 20. 

o nun  f o 


We  note  that  when  R = R . , 3 becomes  3 . The  case  3=0  corresponds  to 

min’  o 

the  lower  limit  of  admissible  R,  R = X /X^. 

’ of 

We  then  get  from  (C3) 

3(3  +X  )F*  (3)  - X F(3)  = 0 (C6) 

o o 


and  from  the  representation  of  (C2) 


- (2F*  (3)  + (3  + X )F"  (3) ) >0 
o 


If  Eq.  (C6)  possesses  no  positive  beta-solution,  then  corresponds 

to  the  lower  limit  of  its  admissible  range,  and  hence  we  have  3^  = 0-  If 

a unique  positive  3-solution  to  Eq.  (C6)  exists,  then  this  is  the  beta  value 

which  determines  and  consequently  3^  is  the  unique  solution*.  If 

Eq.  (C6)  possesses  multiple  solutions  they  must  be  substituted  into  Eq.  (C7) 

and  the  beta  values  which  do  not  violate  (C7)  form  a set  of  admissible  beta 

solutions.  To  this  solution  set  we  must  add  the  lower  limiting  value  3=0. 

We  then  evaluate  Y(X;R)  for  each  of  the  beta  values  of  this  modified  solution 

set.  The  beta  value  which  results  in  the  smallest  Y- value  is  B • 

o 


* 

These  two  cases  apply  for  the  particular  F(X)  functions  presented  in 
Appendix  A. 


1 


For  the  Y-value  comparison  test,  it  is  sufficient  to  consider  any  value 


of  X > and  any  S.  which  satisfies  S.  > 2 0.  By  choosing  S,  = 0 and 

f f if  f 


considering  the  limiting  value  X the  comparison  process  is  greatly 


simplified.  Substituting  beta  into  y(oo,  R)  and  dividing  by  S.X  yields  the 


following  function  of  beta  for  the  comparison  test 


cp(3)  = + \)  F*  (&))  = 


When  multiple  solutions  to  Eq. (C6)  exist,  we  proceed  by  calculating  cp  for 


the  beta  values  of  the  modified  solution  set.  The  beta  which  results  in 


the  smallest  value  of  cp  is 


Once  we  have  obtained  the  value  of  g , Eg. (C5)  allows  us  to  determine 

o' 


R . = (g  + X )/X, 

min  o of 


Ecjuations  (c6 ) and  (C7)  show  a parametric  dependence  on  Should  we 


decide  on  a different  cjuality  factor  Q and  hence  we  must  repeat  the 


process  to  obtain  R . related  to  this  new  \ . 

min  o 


Next  we  consider  for  illustration  the  Y(X)  function 


relation  which  result  from  the  hyperbolic  tangent  function  for  F(X);  hence 


F (g)  = tanh  (g) 


and  (C6 ) becomes 


g(g  + X ) sech  (g)  - X tanh(g)  = 0 
o o 


We  now  select  Q = . 005  and  X = 3 so  that 

o 


2g(g  + 3)  - 3 sinh(2g)  = 0 


The  unicjue  solution  is  g = .4777  and  consequently 


R . = 3.4777/X^;  X = 3 

min  f o 


f 


C4. 


If  we  instead  change  the  quality  so  that  Q = .0009,  then  - 3.853 
and  Eq.  (Cll)  has  the  unique  solution  3 = .3783;  then 


R . = 4.2313/X-;  = 3.853  . 

min  t o 


(C14) 


We  see  from  Eq.  (Cll)  that  for  all  admissible  we  always  have  a unique 
positive  p^.  This  is  similarly  the  case  when  F(X)  is  the  error  function. 


Y(R) 


Fig. 7 
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Chief  of  Naval  Research 
Department  of  the  Navy 
Arl logon,  Virginia  22217 
Attn:  Code  474  (2) 

471 

222 

Director 

ONR  Branch  Office 
495  Summer  Street 
Boston,  Massachusetts  02210 

Director 

ONR  Branch  Office 
219  S.  Dearborn  Street 
Chicago,  Illinois  60604 


Cotrmandlng  Officer 
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Attn:  Mr.  R.  Shea 

U.S.  Anny  Materials  Res.  Agency 

Watertown,  Massachusetts  02172 

Watervliet  Arsenal 
MAGGS  Research  Center 
Watervliet,  New  York  12189 
Attn:  Director  of  Research 

Technical  Library 


Redstone  Scientific  Info.  Center 
Chief,  Document  Section 
U.S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 


Director 

Naval  Research  Laboratory 
Attn:  Code  2629  (ONRL) 

Washington,  D.C.  20390  (6) 

U.S.  Naval  Research  Laboratory 
Attn:  Code  2627 

Washington,  D.C.  20390 

Director 

ONR  - New  York  Area  Office 
715  Broadway  - 5th  Floor 
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Director 

ONR  Branch  Office 
1030  E.  Green  Street 
Pasadena,  California  91101 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  (12) 

Army 


Army  RSD  Center 

Fort  Belvoir,  Virginia  22060 

Commanding  Officer  and  Director 
Naval  Ship  Research  S Development  Center 
Bethesda,  Maryland  20034 
Attn:  Code  042  (Tech.  Lib.  Br. ) 
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(Dr.  W.D.  Sette) 
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(Dr.  M, 

.M.  Sevik) 
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. Strassberq 
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(Dr.  0 Pelt) 

Naval  Weapons  Laboratory 
Dahlgren,  Virginia  22448 


Coirmanding  Officer 

U.S.  Army  Research  Office  Durham 

Attn:  Mr.  J.  J.  Murray 

CRD-AA-IP 

Box  CM,  Duke  Station 

Durham,  North  Carolina  27706  2 


Naval  Research  Laboratory 
Washington,  D.C.  20375 
Attn : Code  8400 

8410 
8430 
8440 
6300 
6390 
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Undersea  Explosion  Research  Div. 
Naval  Ship  R&D  Center 
Norfolk  Naval  Shipyard 
Portsmouth,  Virginia  23709 
Attn:  Dr.  E.  Palmer 

Code  780 


Naval  Ship  Research  S Development  Center 
Annapolis  Division 
Annapolis,  Maryland  21402 
Attn:  Code  2740  - Dr.  Y.F.  Wang 

28  - Mr.  R.J.  Wolfe 
281  - Mr.  R.B.  Niederberger 
2814  - Dr.  H.  Vanderveldt 


Technical  Library 

Naval  Underwater  Weapons  Center 

Pasadena  Annex 

3202  E.  Foothill  Blvd. 

Pasadena,  California  91107 


U.S.  Naval  Weapons  Center 
China  Lake,  California  93557 
Attn:  Code  4062  - Mr.  W.  Werback 

4520  - Mr.  Ken  Bischel 


Commanding  Officer 

U.S.  Naval  Civil  Engr.  Lab. 

Code  L31 

Poi't  Hueneme,  California  93041 


Technical  Director 

U.S.  Naval  Ordnance  Laboratory 

White  Oak 

Silver  Spring.  Maryland  20910 


Technical  Director 

Naval  Undersea  R&D  Center 

San  Diego,  California  92132 


Supervisor  of  Shipbuilding 
U.S.  Navy 

Newport  News,  Virginia  23607 


Technical  Director 
Mare  Island  Naval  Shipyard 
Vallejo,  California  94592 


U.S.  Navy  Underwater  Sound  Ref.  Lab. 
Office  of  Naval  Research 
P.O.  Box  8337 
Orlando,  Florida  32806 


Chief  of  Naval  Operations 
Dept,  of  the  Navy 
Washington,  D.C.  20350 
Attn:  Code  0p07T 


Stra-tegic  Systems  Project  0fri(e 
Department  of  the  Navy 
Washington,  D.C.  20390 
Attn:  NSP-OOl  Chief  Scientist 


Deep  Submergence  Systems 
Naval  Ship  Systems  Command 
Code  39522 

Department  of  the  Navy 
Washington,  D.C.  20360 


Engineering  Dept. 

U.S.  Naval  Academy 
Annapolis,  Maryland  21402 


Naval  Air  Systems  Command 
Dept,  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NAVAIR  5302  Ae^o  & Siructures 

5308  Structures 
52031 F Materials 
604  Tech.  Library 

320B  Structures 
Director,  Aero  Mechanics 
Naval  Air  Development  Center 
Johnsville 

Warminster,  Pennsylvania  189/a 


Technical  Director 

U.S.  Naval  Undersea  R&D  Center 

San  Diego,  California  92132 


Engineering  Department 
U.S.  Naval  Academy 
Annapolis,  Maryland  21402 


Naval  Facilities  Engineering  Cnmnand 
Dept,  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NAVFAC  03  Research  & Dovelo'i- 
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035  Weapons  Dynamics 


Naval  Ship  Engineering  Center 
Prince  George's  Plaza 
Hyattsville,  Maryland  20782 
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Air  Force 
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Attn:  Code  WWRMDD 

AFFDL  (FDDS) 
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WLRC,  Research  Division 
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Attn:  Mechanics  Div. 
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Structures  Research  Division 
National  Aeronautics  & Space  Admin. 
Langley  Research  Center 
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Research  & Technology 
Washington,  D.C.  02546 
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Marine  Corps  Dev.  & Educ.  Command 
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Director 

National  Bureau  of  Standards 
Washington,  D.C.  20234 
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National  Science  Foundation 
Engineering  Division 
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Science  & Tech.  Division 
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Washington,  D.C.  20540 

Director 
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Director  Defense  Research  & Engrg 
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Federal  Avaiation  Agency 
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Maritime  Administration 
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Prof.  R.B.  Testa 
Columbia  University 
Dept,  of  Civil  Engineering 
S.W.  Mudd  Bldg. 

New  York,  N.Y.  10027 

Prof.  H.  H.  Bleich 
Columbia  University 
Dept,  of  Civil  Engineering 
Amsterdam  & 120th  St. 

New  York,  N.Y.  10027 

Prof.  F.L.  DiMaggio 
Columbia  University 
Dept,  of  Civil  Engineering 
616  Mudd  Building 
New  York,  N.Y.  10027 

Prof.  A.M.  Freudenthal 
George  Washington  University 
School  of  Engineering  & 

Applied  Science 
Washington,  D.C.  20006 

D.  C.  Evans 
University  of  Utah 
Cotnputer  Science  Division 
Salt  Lake  City,  Vfesh  Bjjll? 

Prof.  Norman  Jones 
Massachusetts  Inst,  of  Technology 
Dept,  of  Naval  Architecture  S 
Marine  Engrng 

Cambridge,  Massachusetts  02139 

Professor  Albert  I.  King 
Biomechanics  Research  Cent<ir 
lifeyne  State  University 
Detroit,  Michigan  1*9202 


Dr.  V.  R.  Hodgson 
Wbyne  State  University 
School  of  Medicine 
Detroit,  Michigan  1*9202 


Dean  B.  A.  Boley 
Northwestern  University 
Technological  Institute 
2145  Sheridan  Road 
Evanston,  Illinois  60201 


1 


1 


Prof.  P.G.  Hoiige,  Jr. 

Univorsity  of  Minnesota 

Dept,  of  Aerospace  Engng  & Mechanics 

Minneapolis,  Minnesota  55455 

Or.  D.C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbana,  Illinois  61801 

Prof.  N.M.  Newmark 
University  of  Illinois 
I Dept,  of  Civil  Engineering 

I Urbana,  Illinois  61801 

I * 

I Prof.  E.  Reissner 

I University  of  California,  San  Diego 

' Dept,  of  Applied  Mechanics 

La  Jolla,  California  92037 

Prof.  William  A.  Nash 
University  of  Massachusetts 
Dept,  of  Mechanics  & Aerospace  Engng. 
Amherst,  Massachusetts  01002 

I Library  (Code  0384) 

I U.S.  Naval  Postgraduate  School 

Monterey,  California  93940 

j Prof.  Arnold  Allentuch 

Newark  College  of  Engineering 
Dept,  of  Mechanical  Engineering 
323  High  Street 
Newark,  New  Jersey  07102 

Dr.  George  Herrmann 
Stanford  University 
Dept,  of  Applied  Mechanics 
Stanford,  California  94305 

Prof.  J.  D.  Achenbach 
Northwestern  University 
Dept,  of  Civil  Engineering 
Evanston,  Illinois  60201 

Director,  Applied  Research  Lab. 
Pennsylvania  State  University 
P.  0.  Box  30 

State  College,  Pennsylvania  16801 

Prof.  Eugen  J.  Skudrzyk 
Pennsylvania  State  University 
Applied  Research  Laboratory 
Dept,  of  Physics  - P.O.  Box  30 
State  College,  Pennsylvania  16801 


Prof.  J.  Kenipner 

Polytechnic  Institute  nf  Dro(.kl>' 

Dept,  of  Aero.Engrg.E-  /vpplird  Mech 
333  Jay  Street 
Brooklyn,  N.Y.  11201 

Prof.  J.  Klosner 

Polytechnic  Institute  of  Brooklyn 
Dept,  of  Aerospace  & Appl . Merh. 

333  Jay  Street 
Brooklyn,  N.Y.  11201 

Prof.  R.A.  Schapery 
Texas  A&M  University 
Dept,  of  Civil  Engineering 
College  Station,  Texas  77840 

Prof.  W.D.  Pi  1 key 
University  of  Virginia 
Dept,  of  Aerospace  Engineering 
Charlottesville,  Virginia  22903 

Dr.  H.G.  Schaeffer 
University  of  Maryland 
Aerospace  Engineering  Dept. 

College  Park,  Maryland  20742 

Prof.  K.D.  Willmert 
Clarkson  College  of  Technology 
Dept,  of  Mechanical  Engineering 
Potsdam,  N.Y.  13676 

Dr.  J.A.  Stricklin 
Texas  ASM  University 
Aerospace  Engineering  Dept. 

College  Station,  Texas  77843 

Dr.  L.A.  Schmit 

University  of  California,  LA 

School  of  Engineering  S Applied  Science 

Los  Angeles,  California  90024 

Dr.  H.A.  Kamel 
The  University  of  Arizona 
Aerospace  S Mech.  Engineering  Dept. 
Tucson,  Arizona  85721 

Dr.  B.S.  Berger 
University  of  Maryland 
Dept,  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Prof.  G.  R.  Irwin 

Dept,  of  Mechanical  Engrg. 

University  of  Maryland 
College  Park,  Maryland  20742 


P 


L 


r 


Dr.  S.J.  Fenves 
Carnegie^Mellon  University 
Dept,  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Ronald  L.  Huston 

Dept,  of  Engineering  Analysis 

Mail  Box  112 

University  of  Cincinnati 
Cincinnati , Ohio  U5221 


Prof.  George  Sih 
Dept,  of  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 

Prof.  A.S.  Kobayashi 

! University  of  Washington 

Dept,  of  Mechanical  Engineering 
j Seattle,  Washington  98105 

! Librarian 

j Webb  Institute  of  Naval  Architecture 

I Crescent  Beach  Road,  Glen  Cove 

I Long  Island,  New  York  11542 

I Prof.  Daniel  Frederick 

Virginia  Polytechnic  Institute 
Dept,  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 

Prof.  A.C.  Eringen 

Dept,  of  Aerospace  & Mech.  Sciences 

Princeton  University 

; Princeton,  New  Jersey  08540 

Dr.  S.L.  Koh 

School  of  Aero.,  Astro. & Engr.  Sc. 
Purdue  University 
Lafayette,  Indiana  47907 

Prof.  E.H.  Lee 
Div  of  Fngrg.  Mechanics 
Stanford  University 
Stanford,  California  94305 

Prof.  R.D.  Mind! in 
Dept,  of  Civil  Engrg. 

Columbia  University 
S.W.  Mudd  Building 
New  York,  N.Y,  10027 

i 

) _ 


1 


1 

1 

j 

Prof.  S.B.  Dong 

University  of  California  ] 

Dept,  of  Mechanics  ' 

Los  Jingeles  Cal  ifornia  90024  i 

Prof.  Burt  Paul  • 

University  of  Pennsylvania  ; 

Towne  School  of  Civil  & Mech.  Engrg.  j 

Rm.  113  - Towne  Building  | 

220  S.  33rd  Street 

Philadelphia,  Pennsylvania  19104  j 

Prof.  H.W.  Liu  | 

Dept,  of  Chemical  Engr.  % Metal. 

Syracuse  University 
Syracuse,  N.Y.  13210 

Prof.  S.  Bodner 

Technion  R&D  Foundation  ' 

HajjTi,  Israel 

Prof.  R.J.H.  Bollard 

Chairman,  Aeronautical  Engr.  Dept. 

207  Guggenheim  Hall 

University  of  Washington  j 

Seattle,  Washington  98105 


Prof.  G.S.  Heller 
Division  of  Engineering 
Brown  University 
Providence,  Rhode  Island  02912 

I 

i 

Prof.  Werner  Goldsmith  | 

Dept,  of  Mechanical  Engineering  i 

Div.  of  Applied  Mechanics 
University  of  California 
Berkeley,  California  94720 

Prof.  J.R.Rice 
Division  of  Engineering 
Brown  University 
Providence,  Rhode  Island  02912 

Prof.  R.S.  Rivlin 

Center  for  the  Application  of  Mathematics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 

Library  (Code  0384) 

U.S.  Naval  Postgraduate  School 
Monterey,  California  93940 

« 

Dr.  Francis  Cozza^elli  • 

Div.  of  Interdisciplinary 
Studies  & Research 

School  of  Engineering  f 

State  University  of  New  York 
Buffalo,  N.Y.  14214 


Industry  and  Research  Institutes 


FT 


Library  Services  Department 
Report  Section  Bldg.  14-14 
Argonne  National  Laboratory 
9700  S.  Cass  Avenue 
Argonne,  Illinois  60440 

Or.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
I 129  Mount  Auburn  St. 

I Cambridge,  Massachusetts  021 3B 

I 

I Dr.  L.H.  Chen 

' ■ General  Dynamics  Corporation 

I Electric  Boat  Division 

Groton,  Connecticut  06340 

Dr.  J.E.  Greenspon 

J.G.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 

. Dr.  S.  Batdorf 

I The  Aerospace  Corp. 

P.O.  Box  92957 

Los  Angeles,  California  90009 

! Dr.  K,C.  Park 

Lockheed  Palo  Alto  Research  Laboratory 
Dept.  5233,  Bldg.  205 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

Library 

Newport  News  Shipbuilding  & 

Dry  Dock  Company 

Newport  News,  Virginia  23607 

, Dr.  W.F.  Bozich 

McDonnell  Douglas  Corporation 
5301  Bolsa  Ave. 

Huntington  Beach,  CA  92647 

Dr.  H.N.  Abramson 

Southwest  Research  Institute 

Technical  Vice  President 

Mechanical  Sciences 

P.O.  Drawer  28510 

San  Antonio,  Texas  78284 


1 

=* 


Or.  R.C.  DeHart 

Southwest  Research  Institute 

Dept,  of  Structural  Research 

P.O.  Drawer  28510 

San  Antonio,  Texas  78284 

Dr.  M.L.  Baron 
Weidlinger  Associates, 

Consulting  Engineers 
110  East  59th  Street 
New  York,  N.Y.  10022 

Dr.  W.A.  von  Riesemann 
Sandia  Laboratories 
Sandia  Base 

Albuquerque,  Ney  Mexico  87115 
Dr.  T.L.  Geers 

Lockheed  Missiles  8 Space  Co. 

Palo  Alto  Research  Laboratory 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Dr.  J.L.  Tocher 

Boeing  Computer  Services,  Inc. 

P.O.  Box  24346 

Seattle,  Washington  98124 

Mr.  Wi Ilian  Caywood 

Code  BBE,  Applied  Physics  Laboratory 

8621  Georgia  Avenue 

Silver  Spring,  Maryland  20034 

Mr.  P.C.  Durup 
Lockheed-California  Company 
Aeromechanics  Dept.,  74-43 
Burbank,  California  91503 


Unclassified 

security  classification  of  this  page  f»Ti»n  Dmim  EntaredJ 

REPORT  DOCUMENTATION  PAGE  befoII'^com^le^S^form 

1.  RCR^OT  NUMBSa — - l2.  GOVT  ACCESSION  NO.  3 RECIPIENT’S  CATALOG  NUMBER 


I.  RCPORT  MUMBBa — 


/^(  RPI-CS-77-1 


,4.  TJTLE.CAOtf  Sulx((/Jj  5.>.»VPe  OF  REPORT  « PER4O0  COVERED 

Nonlinear  Monotonic  Functions  with  Selectable  Topical  Report*^ 

Intervals  of  Almost  Constant  or  Linear  Behavior  ~ ' 

with  Application  to  Total  Strain  Viscoplasticity*  ® performing  org.  report  number 


Topical  Report*^ 


17.  AUTHORf»J 


E.P.  Cernocky  and  E.  Krempl 


19.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  & Mechanics 
Rensselaer  Polytechnic  Institute.  Troy.  NY  12181 


».  CONTRACT  OR  5,SAaI  MUMBERf»> 

NC^0014-76-C-0231  ( ^ 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


NR  064-571 

12.  PEPORT  OATe” 


II.  CONTROLLING  OFFICE  NAME  and  ADDRESS  12  REPORT  DATE 

Dept,  of  the  Navy,  Office  of  Naval  Research  / ) januaav  1977 

Structural  Mechanics  Program  is.  nu»^r 

Arlington,  Va.  22217  /,  / J.  " , 

14.  MONITORING  AGENCY  NAME  A ADDRESS^ff  dlllarani  from  ControtUnt  Otilct)  15.  SECURITY  CljASS.  (ol  (h(«  reportj 

Office  of  Naval  Research  - Resident  Representative  Unclassified 

715  Broadway  - 5 th  Floor  ^ 

New  York,  NY  10003  declassification  downgrading 


16.  DISTRIBUTION  ST  ATEMENT  fof  (M«  ReporlJ 


Approved  for  public  release;  distribution  unlimited. 


I 17.  DISTRIBUTION  STATEMENT  (oi  the  abatract  antarad  In  Block  20,  It  dlHarant  from  Raport) 


I 18.  SUPPLEMENTARY  NOTES 


19  KEY  WORDS  (Continua  on  ravaraa  aide  It  nacaaaary  and  Identify  by  block  number) 

Plasticity,  viscoplasticity,  cyclic  loading,  rate  effects,  analytic 
modelling,  numerical  experiments,  curve  fitting  of  mechanical  behavior 
data. 


I 20.  ABSTRACT  ^Confinuo  on  ravaraa  aide  If  nacaaaary  and  Identify  by  ftfocit  number) 


“pR*  A rather  general  method  is  given  to  construct  classes  of  functions  with 
an  arbitrary  almost  constant  (linear)  initial  interval  followed  by  a non- 
prescribed  interval  of  monotonic  nonlinear  behavior.  This  region  of  nonlinear 
behavior  is  succeeded  by  an  unbounded  interval  of  almost  constant  (linear) 
behavior.  They  contain  not  more  than  four  selectable  parameters  and  are 
synthesized  from  analytic,  monotonic,  normalized  and  bounded  base  functions 
through  the  introduction  of  two  separate  kernel  sets,  subsequent  addition  and 


DD  ,:?n7,  1473  EDITION  OF  I NOV  65  1$  OBSOLETE 

S/N  OtO}-OI4-  6601  I 


Unclassified  / 

security  CLASSIFICATION  OF  TRIS  RAGE  fBTi»n  Data  BnttH) 


sif ied 


)W|Ty  classification  JF  THIS  PAOErKhen  Data  EntareJ) 


> integration.  As  examples  we  give  the  special  functions  based  on  the  error, 
the  hyperbolic  tangent,  the  inverse  tangent,  a rational  and  the  incomplete 
gamma  function.  Limiting  function  forms,  such  as  the  bilinear  form,  are 
derived  for  limiting  values  of  the  parameters. 

We  have  found  these  functions  useful  in  the  total  strain  approach  to 
viscoplasticity,  i.e.,  the  analytical  modelling  of  stress-strain  diagrams, 
strain  (stress) -rate  effects,  creep  and  relaxation  curves  for  monotonic  and 
cyclic  loading.  Also  these  functions  offer  great  flexibility  in  the  curve 
fitting  of  experimental  data  generated  in  the  above-mentioned  tests. 


Unclassified 


